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Abstract. We discuss the role of feedback via photoionization and Compton heating 
in the co-evolution of massive black holes at the center of spheroidal galaxies and their 
stellar and gaseous components. We first assess the energetics of the radiative feedback 
from a typical quasar on the ambient interstellar gas. We then demonstrate that the 
observed Mbh-ct relation could be established at a relatively early epoch in galactic 
evolution when the formation of the stellar bulge was almost completed and the gas- 
to-stars mass ratio was reduced to a low level ~0.01 such that cooling could not keep 
up with radiative heating. A considerable amount of gas was expelled at that time and 
black hole accretion proceeded at a much lower rate thereafter. 

1 General Picture 

Most elliptical galaxies are poor with respect to interstellar gas. Also, elliptical 
galaxies invariably contain central massive black holes (BHs), and there exists 
a tight relationship between the characteristic stellar velocity dispersion a and 
the BH mass Mbh [1,2], and between A/bh and the host spheroid mass in stars, 
Af, [3]. Are these two facts related? Here we focus on a scenario in which the 
mass of the central BH grows within gas rich elliptical progenitors until star 
formation has reduced the gas fraction in the central regions to of order 1% of the 
stellar mass. Then radiative feedback during episodes when the BH luminosity 
approaches its Eddington limit drives much of the gas out of the galaxy, limiting 
both future growth of the BH and future star formation to low levels. 

Many works already recognized the importance of feedback as a key ingredi- 
ent of the mutual BH and galaxy evolution [4-16]. What is new about this work 
is the stress on one component of the problem that has had relatively little at- 
tention: the radiative output of the central BH is not conjectural - it must have 
occured - and the high energy component of that radiative output will have a 
dramatic and calculable effect in heating the gas in ellipticals. 

Using the average quasar spectral output derived in [17], we show below and 
in more detail in [18] that the limit on the central BH produced by the above 
argument coincides accurately with the observed A/bh""" relation. Not only the 
slope, but also the small normalization factor is nicely reproduced. 
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The present work is complementary to [5,8] in that, while it does not attempt 
to model the complex flaring behavior of an accreting BH with an efficient hydro- 
dynamical code, it does do a far more precise job of specifying the input spectrum 
and the detailed atomic physics required to understand the interaction between 
that spectrum and the ambient interstellar gas in elliptical galaxies. 

2 Radiative Heating of ISM in Spheroids 

Below we assess the conditions required for the central BH radiation to signifi- 
cantly heat the ISM over a substantial volume of the galaxy. In this section we 
shall assume that the central BH has a mass as given by the observed Mbh-c 
relation for local ellipticals and bulges [2]: 

MBH = 1.5xlO«Me(^^^^)\ (1) 

This assumption will be dropped in §3, where we predict the Mbr-ct relation. 

In [17] we computed for the average quasar spectrum the equilibrium tem- 
perature Teq (at which heating due to Compton scattering and photoionization 
balances cooling due to line and continuum emission) of gas of cosmic chemical 
composition as a function of the ionization parameter ^ = L/nr^, where L is the 
BH bolometric luminosity. In the range 3 x 

Teq(O^2xl02^K, (2) 

while at C < 10^ and ^ > 10^ Teq « 10'' K and 2 x lO'^K, respectively. On the 
other hand, the galactic virial temperature is given by 

^--^•0^^«'^(2001^)'- 
We can then find the critical density ncrit defined by 

Teq(i/ncritr^) = Tvir (4) 

as a function of distance r from the BH. Gas with n < ncrit('') will be heated 
above Tvir and expelled from the galaxy. We show in Fig. 1 the resulting (r, n) 
diagrams for a small and large BH/galaxy. 

In reality, provided that Tgq > Tvir, significant heating will take place only 
out to a certain distance that depends on the luminosity and duration of the 
quasar outburst. Since the BH releases via accretion a finite total amount of 
energy, cMbhc^, there is a characteristic limiting distance: 

Inside this radius, a low density, fully photoionized gas will be heated to the 
Compton temperature Tq « 2keV characteristic of the quasar spectral output. 
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Fig. 1. The (r, n) plane for a galaxy with a = 180 kms"^ (Tvir = 2.4 x 10® K, Mbh = 

10* Mq, upper panel), and with a = 320 kms^^ (Tvir = 7.7 x 10® K, A/bh = IO^Mq, 
lower panel). In the dashed area, gas can be heated above Tvir by radiation from the 
central BH emitting at the Eddington luminosity. The upper boundary of this area 
scales linearly with luminosity. Vertical boundaries are Rc, Ri, R2 and R^. 



More relevant for the problem at hand is the distance out to which low density 
gas will be Compton heated to T ^ 

=l,300pc(-^) — — ^ (6) 

Tvir/ VO.I/ 200 knis-i ^ ' 

Yet another characteristic radius is out to which gas of critical density ricrit will 
be heated to T ^ Tyir via photoinization and Compton scattering: 

R2 = Ri[r{n„i^)/rc]'/\ (7) 

where /c and F are the Compton and total heating rates, respectively. Depend- 
ing on gas density (0 < n < Ucrit), the outer boundary of the "blowout region" 
will be located somewhere between Ri and i?2- The size of the heating zone can 
be compared with the galaxy effective radius 

^^-^'000P<2001^)'- 

The different characteristic distances defined above are shown as a function 
of Mbh in Fig. 2. One can see that a BH of mass < IO^Mq should be able to 
unbind the ISM out to several R^. In the case of more massive BHs/galaxies with 
Mbh ~ IO^-IO^Mq, the heating will be localized to innermost ~0.3-0.5iie- 
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Fig. 2. Different heating radii: Rc (dotted line), Ri (short-dashed hne), and R2 (long- 
dashed line), and the galactic effective radius (solid line), as a function of Mbh- 

3 Possible Origin of the Mbh— cr Relation 

We now consider the following general idea. Before the BH grows to a certain 
critical mass, MBH,crit, its radiation will be unable to efficiently heat the ambient 
gas, and the BH will accrete gas efficiently. Once the BH has grown to MBH.crit; 
its radiation will heat and expell a substantial amount of gas from the central 
regions. Feeding of the BH will then become sclf-rcgulatcd on a time scale of 
order the cooling time of the low density gas. Subsequent central activity will 
be characterized by a very small duty cycle (~0.001), as predicted by hydrody- 
namical simulations [5,8] and suggested by observations [14]. BH growth will be 
essentially terminated. 

Suppose that the galaxy density distribution is that of a singular isothermal 
sphere, with the gas density following the total density: 



Here Mgas and M are the gas mass and total mass within the region affected by 
radiative heating. The size of the latter is uncertain but is less than a few kpc 
(see §2), so that M is dominated by stars rather than by dark matter. 
Radiation from the central BH can heat the ambient gas up to 



this approximate relation being valid in the range 3 x 10^-10^ K. Remarkably, 
Teq does not depend on distance for the adopted density distribution. We 




(9) 




(10) 
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then associate the transition from rapid BH growth to slow, feedback-limited 
BH growth with the critical condition 

-^eq ^esc-^ir? (-^-^) 

where r^csc <; 1 and Tvir is given by (3). Once heated to Teq ^ Tvir, the gas will 
stop feeding the BH. The condition (11) will be met for 

Me„„,. = 4.6 X IO.»M,,.„ (_£_,)* iSl^fc. 

Therefore, for fixed values of r/eso L/L-^dd and Mg^/M we expect MBH,crit cx a'*, 
similarly to the observed Mbu-ct relation. Equally important is the normaliza- 
tion of the A^BH o" relation. By comparing (12) with (1) we find that the observed 
correlation will be established if 

Mgas/M = 3 X 10-3,7- jL/LEdd- (13) 

The gas-to-stars ratio is thus required to be low and approximately constant for 
spheroids of different mass at a certain stage of their evolution. As for the Ed- 
dington ratio, it is reasonable to expect L/L^dd ~0.1-1 during quasar outbursts. 




tJg (km/s) 

Fig. 3. Thick solid line shows the predicted Mbh-c correlation resulting from heating 

of the ISM by the radiation from the central BH assuming Mg^s{Rc)/M — 0.003 and 
»7esc = 1. Thin solid line corresponds to Mgas{Re)/M = 0.0015 and ?7osc — 2. Dashed 
line is the observed Mbh oc a'^ correlation in the range 10^-a few 10® M©, extrapolated 
to lower and higher BH masses. Dotted lines are Mbh oc and Mbh oc laws. 

The approximately linear Teq(^) dependence [see (2)] was crucial to the above 
argument leading to the MeH.crit oc cr^ result. However, the Teq(^) function 
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becomes nonlinear outside the range 3 x lO^K < T^q < lO'^K [17]. In Fig. 3 

we show the predicted correlation between A^BH.crit and a for L/L-^dd = 1 and 
Mgi^s/M = 3 X 10^^. It can be seen that the Mbh oc ct'* behavior is expected 
to break down for Mbh < W^Mq and for Mbh ^ a few 10® Mq. It is perhaps 
interesting that the range of masses shown in Fig. 3 for which A^bh k cr^ is 
obtained from considerations of atomic phyisics (and the observed AGN spectra) 
corresponds closely with the range of masses for which this power law provides a 
good fit to the observations. Exploring the A^bh o" relation observationally near 
IO^Mq would be a sensitive test of the importance of radiative feedback. 

4 Detailed Modelling of the BH-Galeixy Co-evolution 

In [18-20] we addressed in a more quantitative way the BH growth in the context 
of the parent galaxy evolution. We adopted a physically-motivated one-zone 
model, taking into account the mass and energy return from the evolving stellar 
population. This model predicts that after an initial "cold" phase dominated by 
gas infall, once the gas density becomes sufficiently low the gas heating dominates 
and the galaxy switches to a "hot" solution. The gas mass/stellar mass ratio at 
that epoch (~0.003) is remarkably close to the value inferred above from the 
argument leading to the right Mbh relation. Other predictions of the toy 
model are also in satisfactory agreement with observations. The "cold" phase 
would probably be identified observationally with the Lyman Break and SCUBA 
galaxies, while the "hot" phase with normal, local ellipticals. 

A proper investigation of the importance of radiative heating on the BH/galaxy 
co-evolution, based on hydrodynamical numerical simulations, is now in progress. 
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